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Abstract 

We consider lattice gauge theories at strong coupling with gauge group U(Nc), or SU(Nc) 
restricted to the meson sector, and coupled to Nf flavors of fundamental representation 
staggered fermions. We study the formation of a chiral condensate by means of resummation 
of a hopping expansion. Different classes of graphs become dominant as the parameter 
(Np/Nc) is varied. By performing graph resummation we obtain an equation for determining 
the condensate as a function of (Np/Nc) and mass m. For values of {Np/Nc) below a critical 
value one reproduces the well-known result of the existence of a non-vanishing condensate 
solution in the m = limit. Above the critical {Np/Nc) value, however, no such solution 
exists, its abrupt disappearance indicating a first order transition to a chirally symmetric 
phase with composite (colorless) excitation spectrum. 
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1 Introduction 



Chiral symmetry breaking in QCD is one of the cornerstones of modern particle physics. 
The existence of a non-vanishing chiral condensate at T = has of course been established 
and studied in any number of Monte Carlo investigations over several decades pQ. The 
phenomenon can be nicely demonstrated analytically in lattice QCD at strong coupling as 
was done for large number of colors Nc long ago [2] - [3] . 

Conventional wisdom seems to have been that chiral symmetry will eventually be broken 
by taking the coupling sufficiently strong regardless of the number of fermion flavors. The 
reason one might suspect this is not correct is discernible in the very same analytic demon- 
strations of chiral symmetry breaking at strong coupling: as the number of flavors becomes 
large compared to the number of colors, contributions in the fermion hopping expansion that 
used to be subdominant grow to be dominant, and the original result of a non-vanishing 
chiral condensate is no longer necessarily valid. 

In this paper we use the same framework of lattice gauge theory at strong coupling and 
gauge group U{Nq), or SU(Nc) restricted to the mesonic sector, but with any number Np 
of flavors of fundamental representation staggered fermions, and employ a hopping expansion 
to study the formation of a chiral condensate. To see whether a condensate forms in the 
chiral limit m — > requires resummation of infinite classes of graphs allowing one to obtain 
an equation for the condensate. One then finds that the outcome depends on the parameter 
(Np/Nc)- There is a critical value of this parameter below which the usual non- vanishing 
condensate forms. Above this critical value, however, the non-vanishing condensate solution 
abruptly ceases to exists as a physical solution (turns complex), which is the usual sign of a 
first order phase transition. 

Such a result has in fact already been demonstrated for the case Nc = 3 in a recent 
Monte Carlo study [5J. These authors also provide evidence that the chiral phase resulting 
above the critical number of flavors possesses a massless spectrum. 

It should be noted that it is perfectly straightforward to set up a 1/Np expansion in the 
continuum theory by, as usual, keeping a 't Hooft coupling A 2 = Npg 2 fixed and appropriately 
rescaling the gauge field and fermion sources. Such an expansion [6], however, is not useful for 
the problem addressed here, i.e., the strong coupling regime, as it is incapable of providing 
even a hint of the actual spectrum. For that one has to use the nonperturbative lattice 
formulation, to which this continuum expansion set up cannot be carried over, and a rather 
more involved procedure is needed to deal with the Np versus Nc dependence. 

Existence of strongly coupled gauge systems with composite (colorless) spectrum and 
chiral symmetry is of great physical and theoretical interest on many fronts. Apart from 
the basic questions of the complete phase diagram in (/3, (Np/Nc),m), and the relation to 
possible continuum limits, there is a host of potential applications for such systems both 
in particle and condensed matter physics (cf remarks in section [5|). A lot of work will be 
required to adequately explore these various directions. 

The paper is organized as follows. In section [2] we set up the expansion and classify graphs 
according to their dependence in Np, Nc and (inverse) powers in m . This allows us to see 
how to group graphs in classes to be resummed so that the small m limit may eventually be 
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considered. The resummation is then done in section [3] and the equation for the condensate 
obtained. Its solutions are considered in section HI Section [5] contains our concluding remarks. 



2 Hopping expansion and (Np/Nc) graphology 

We work on a euclidean hypercubic d-dimensional lattice with lattice sites denoted by their 
lattice coordinates x = (x 11 ), and lattice unit vectors in the /i-th. direction by ft. We use 
standard lattice gauge theory notations and conventions. The gauge field bond variable £/& 
on bond b = (x,jl) is more explicitly denoted by U^x), and the fermion fields on site x by 
V>(x) and ip(x). 

The lattice action for massless fermions is given by: 

S F = \ [^{x) 1 ^x)U^{x)^{x + (i)-^{x + fi) 1 ^x)U} 1 {x)ij{x)Y (2.1) 

b=(x,fi) 

The matrices 7 [6] = 7^(x) defined on each bond b = (x,fi) satisfy n&edp7M = ^ or eacn 
plaquette p. For staggered fermions 

7m (x) = (-1)^<m^ ) 7 i(x) = 1. (2.2) 

For naive fermions Ja(x) = 7„ are elements of the Euclidean Dirac-Clifford algebra satis- 
fying {7^,7^} = 25^1 and are hermitian matrices. Naive fermions are of course unitarily 
equivalent to 2 d l 2 copies of staggered fermions. 

We take the fermions in Np flavors or tastes. For staggered fermions this corresponds 
to ANp continuum flavors. The staggered fermion action in (|2.ip is then invariant under 
a U{Np) x U{Np) global symmetry corresponding to independent rotations of fermions on 
even and odd sublattices [7j. This symmetry is referred to as chiral symmetry. 

We take the fermions to transform under the fundamental representation of the gauge 
group which is taken to be U(Nc) or SU(Nc)- In the case of SU(Nc), however, we restrict 
to the meson sector in the hopping expansion below. This is not essential and is done to 
avoid the technical nuisance of having to also include the baryon sector which would made 
things unnecessarily more cumbersome. 

We consider the order parameter for chiral symmetry breaking ip(x)'ip(x). Its vev is related 

to the fermion 2-point function (full propag ator) G ai ' bj (x,y) = (^ ai {x)^ {y)^j in the limit 
x = y: 

($(x)ip(x)\ = -tr[G(x,x)} . (2.3) 

Here tr denotes trace over color, flavor, and (if naive fermions are considered) Dirac spinor 
indices. We will denote by tr c and tr F traces over color (a) and flavor (i) indices, respectively. 
Note that G{x, x) = tr c G(x, x) is a gauge invariant quantity. 

To study (|2.3p we add an external source coupled to the operator tp(x)tp(x), i.e. a mass 
term, to the action. We write the lattice gauge theory action in the presence of the source 
more concisely in the form 

s=s G +s F = Y,Pi 1 -^ RetTU p]+J2^ x ) K ^ u ^^ ' ( 2 - 4 ) 
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where 
with 



K x>y {U) = M x>y (U)+K x>y (2.5) 



M x,yi U ) = 2 ^ U i^( x ) S y,x+p,- 1/j. u I( x ~ fr) S y,x-p. (2-6) 
K XiJ , = ml c l F 5 X: y . (2.7) 



In the strong coupling limit /3 — > the plaquette term in (|2.4p is dropped. The corrections 
due to this term can be taken systematically into account within the strong coupling cluster 
expansion, which, for sufficiently small /3, converges. 

Setting (3 = in (|2.4p then, G(x, x), in the presence of arbitrary source K, is given by 

G{*,x) = j^^ flDin^in^ (2.8) 

J[DU] Det[l + K^MjU)] [[1 + K^MjU)]^^ 1 ]^ 

f[DU] Det[l + K- 1 M(C/)] ' ( ' ' 

The vev of ip(x)'ip(x) in the presence of the source is then obtained from (|2.3[) . 



We evaluate (|2.9p in the hopping expansion. This amounts to expanding (|2.9p treating 
M as the interaction and taking K as defining the inverse bare propagator: 

K-J = m- 1 l c l F S x>y . (2.10) 

Note that K is purely local, whereas M has only nearest-neighbor non-vanishing elements 
^-x,x+p, — It^C^u(^) an d Mj._j._j; = —^j^U^x — fj,). The expansion can be shown to converge 
in the large volume limit for sufficiently large m. In the absence of the plaquette term, and 
since M is linear in the bond variables U^, integration over the gauge field results in non- 
vanishing contributions only if at least two M factors with equal number of _Vs and U^'s 
occur on each bond b. In the case of SU(Nc) equal mod Nq number of Ub and U b s are also 
allowed but will be excluded by restriction to the mesonic sector as already explained above. 



x x 
(a) 



(b) 



Figure 1: (a) Some tree graphs; and (b) some loop graphs attached to site x. (Arrows indicate 
U, W'a on bonds.) 



The expansion of the K. XX (U) is represented by all paths starting and ending at x, whereas 
that of the Det/C(f/) by all closed paths [l], [HI- Consistent with the above constraint on each 
bond resulting from the [/-integrations, after the cancellation of all disconnected graphs 
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between numerator and denominator the connected graphs giving the expectation (|2.9[) nat- 
urally fall into two classes: 'tree graphs' and 'loop graphs'. 

The tree graphs consist of paths starting and ending at x and enclosing zero area (Fig. 
[T] (a)); they arise entirely from the expansion of K^\(U). The loop graphs, such as those 
depicted in Fig. QJb), consist of paths from the expansion of 1C~],.(U) and loops from that of 
Det/C(C7) coupled by the [/-integrations in the numerator in (12.9|) [9]. Though a well-known 
fact (see e.g. [T]), it is perhaps worth recalling here that within this hopping expansion 
there are no restrictions on how many times a bond is revisited in drawing all such possible 
connected graphs [ID] . 

This hopping expansion is naturally organized as an expansion in powers of 1/Nq, for 
fixed Np. The trees graphs are all or order unity, whereas loops graphs contribute inverse 
powers of Nc- Note that, given this expansion of G(x,x), the vev of iptp contains an extra 
overall multiplicative factor of NcNp from the trace in (12.31) . 

The trees then give the leading contribution in the large Nq limit. Self-consistent sum- 
mation of all trees leads to the well known result [2] - [3] , [II] of non- vanishing vev fl2.3f> 
in the limit of vanishing external source (m — > 0), i.e. chiral symmetry breaking at strong 
coupling. 

The loop contributions provide non- leading small corrections for large Nc- They, however, 
also contain Np dependence. This can lead to drastically different state of affairs if Np is no 
longer taken to be small compared to Nc- Thus, among the graphs in Fig. [I] (b) the first 
is of order (Np/Nc)-, the second is of order (Np/Nc) 2 , and the third is of order (Np/Nc) 4 ■ 
They are indeed subdominant for large enough Nq at fixed Np; but they become dominant 
in the opposite limit of large (Np/Nc)- To consider this latter limit then we have to deal 
with the loop graphs. 



xxx x 

(a) (b) (c) 

Figure 2: (a) A maximally tiled graph of area 4 attached at site x; (b) two non-maximal 
tilings of the graph; (c) a maximally tiled 2-plaquette graph with one plaquette multiply 
tiled. (Plaquettes in (c) drawn to larger scale for clarity. Also, for simplicity, here and in 
subsequent figures we no longer always indicate the U and directions along each graph 
bond.) 



Consider the (Np/Nc) order of graphs of a given non-zero lattice area A (in plaquette 
units). The graphs of the highest possible order, (Np/Nc) A , are those maximally tiled, i.e., 
tiled by a number of single plaquette loops equal to the area enclosed by the graph. Thus, 
the graph in Fig. EUa) is or order (Np/Nc) 4 ', whereas in Fig. 0(b) the first graph is of order 
(Np/Nc) 3 and the second of order (Np/Nc) 2 - This rule holds whether a graph is planar or 
not (self- intersecting) . If a plaquette is covered k times by the graph folding onto itself, it 
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contributes k units to the area. The same applies in the case of multiply tiled plaquettes; a 
/c-times tiled plaquette is also counted as contributing to the graph area k units [12]. So, e.g., 
in Fig. [2](c) the connected contribution from the 3-tiled plaquette results in order (Np / Nq) 
for the graph, which is then counted as being of area 4. 

We now classify the connected maximally tiled graphs as follows. Unreduced graphs are 
defined as those consisting of a set of tiled plaquettes where each plaquette is connected 
to at least one other plaquette in the set at a common boundary site or along a common 
boundary bond, or is attached to site x. Furthermore, on any plaquette or at site x, tree 
segments of any length and not connecting to any other plaquette may be attached. Thus 
all three graphs in Fig. Q] (b) belong to this class. Reduced graphs are those where each 
tiled plaquette in the graph is connected to one or more of the other plaquettes in the graph 
or to the site x only by a tree segment, with any two plaquettes connected by at most one 
tree segment. Furthermore, on any plaquette or at site x, tree segments of any length and 
not connecting to any other plaquette may be attached. Thus, in common physics jargon, 
reduced graphs are 1PR with respect to any tree segment [13] . Partially reduced graphs then 
are those containing one or more unreduced plaquette subsets connected to the rest of the 
graph or the site x by tree segments. Examples are shown in Fig. [3[ 























X 
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(a) 



(b) 



(c) 



(d) 



(e) 



(f) 



Figure 3: An unreduced graph (a), and one reduced graph (b) and one partially reduced 
graph (c) arising from it; similarly for the unreduced graph (d) and reduced (e) and partially 
reduced (f) graphs. 



Consider all unreduced graphs attached to site x of given total area A and having k 
tree segments of total length I attached (cf. Fig. [3] (a), (d)). Any unreduced graph can be 
constructed by attaching one plaquette at a time either to the site x or to another plaquette 
at a site or along a bond on its boundary. Now, counting all possible orientations, there are 
2d(d— 1) ways of attaching one tiled plaquette at site x. Each subsequent tiled plaquette can 
then be attached in 16(d— 1) ways along a bond, and less than 32d(d— 1) ways at a common 
boundary site, counting all possible positions and flips in every 2-plane. There are 1 + 8 ways 
of attaching the first unit of each segment on a plaquette or at x, i.e, 9 k ways for k segments; 
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and since 1 < k < I, as each segment must be at least one bond long, the total number of 
ways of attaching total-length-/ tree segments is less than 2x9'. Furthermore, there are 
(2d — 1) ways of orienting each segment unit with a backtracking constraint (see below - this 
though does not really affect the argument). So, for the total number of unreduced graphs 
of area A with tree segment of total length I one has 

n u (A, I) < Ad(d - l)[16(d - l)(2d + 1)] A ~ V(2d - 1)' . (2.11) 

Consider now the class of reduced graphs obtained by separating plaquettes by inserting 
tree segments (cf. Fig. [3] (b), (e)). A tiled plaquette can be attached to a tree segment 
endpoint in 8d(d — 1) different ways. Each unit of a connecting segment can be oriented in 
(2d — 1) ways, so if the connecting segments have total length l c , there are (2d — l) lc ways of 
attaching them. Note that there must be A such connecting segments between the plaquettes 
and x, so that l c > A. In addition, k further tree segments of total length I can be attached 
to the A plaquettes or the site x in (1 + 8) k ways. Again, since 1 < k < I, the number of 
ways of attaching total- length-/ tree segments is greater than 9', and they can be oriented in 
(2d — l) 1 ways. Actually, the first leg of segments attached directly at x have no backtracking 
constraint (see below). So the total number of reduced graphs of area A, with connecting 
segments total length l c and other tree segments of total length I attached is not less than 
(d> 3): 

n r (A,l c ,l) > [8d(d-l)] A (2d-l) lc 9 l (2d-l) 1 (2.12) 

> 8d(d - l)(2d - l)[8d(d - l)(2d - l)] A ~ 1 9 l (2d - l) 1 

> n u (A,l) . 

By the same counting, for th.6 number Tip r (A, q, l c , I) of partially reduced graphs with q reduced 
plaquettes, 1 < q < A, and connecting segments of total length (l c — A + q), one has 

n u (A, I) < n pr (A, q, l c , I) < n r (A, l c , , /) . (2.13) 

The reduced graphs then dominate in number among maximally tiled graphs. They also 
clearly contain the highest powers of 1/m. The explicit evaluation of reduced graphs within 
the hopping expansion is straightforward and shows that each reduced tiled plaquette may 
be taken to contribute a factor of (— l)(A^/iVc)(l/4m 2 ) 4 , each tree segment unit (one bond 
long, double line) a factor of (— l/4m 2 ), and the attachment to site x a final factor of (1/m). 
The value J 7 of a reduced graph 7 of given A, l c and I is then independent of how these 
building blocks are arranged, i.e. of the actual reduced graph structure: 

/»r \ A 1 / 1 \ 4A+l c +l 



h = ( -D— ~{^) =U. (2.14) 

Unreduced (q = 0) graphs miss at least a factor of (— l/4m 2 )' c ; this is for unreduced plaquettes 
attached at common boundary sites - unreduced plaquettes attached at common boundary 
bonds miss additional factors of (l/2m) 2 . Partially reduced graphs miss at least a factor 
(— l/4m 2 ) A ~ q . The value I^(A,q,l c ,l) of any unreduced or partially reduced graph 7 can 
then be expressed in terms of the fully reduced graph value (|2.14p as 

I = Cl (-4m 2 ) A - q (4m 2 )^I r , (2.15) 
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where c 7 is a 7-dependent positive or negative numerical constant of order unity, and p 7 > 
a 7-dependent non-negative integer power. Let T r (A,l c ,l) and F pr (A.q,l c ,l) denote the sets 
of reduced and of partially reduced or unreduced graphs, respectively, of given A, l c , I. Then 
by (|2.15p one has 



A-l 



q=0 -y£T P r(A,q,l c ,l) 



,-y£V r (A,l c ,l) 



with 



A-l 



K(A, l c , I) = [n r (A, l c , l)]' 1 E c 7 (-4m 2 ) A - 1 " <? (4m 2 r 

q=0 y£r pr (A,q,l c ,l) 

For m < 1/2 and any constant c > c 7 , use of (]2.13p gives the bound: 



\K(A,l c ,l)\ < c 



1 - (4m 2 ) 



2\A 



(2.16) 



(2.17) 



(2.18) 



1 - (4m) 2 

Consider then the contribution of all maximally tiled graphs up to a maximum size A < K, 



l c <K,l<K: 



K 



K A-l 



E E h + EE E h 

K A,l c ,l 7 eT r (A,l c ,l) ) \A,l c ,l q=0 yeT pr (A,q,l c ,l) 

E E I-y)[l + HK,rn) 

A,l c ,l yeT r (A,l c ,l) J 



(2.19) 



Using (|2.14|) and (|2.16|) one now obtains for the ratio of the contribution of all partially 
reduced graphs to that of all reduced graphs: 



K A-l 



K 



mm)= EE E M E E h 



A,l c ,l g=0 7er pr (A,q,Z c ,Z) 



A,i c ,i^£r r {A,i c ,i) 



M{K) 
V{K) 



where 



K-l 



v = i+ e 



k— (fe-j^ , k<2 .k 3 ) 
k 1 >k 2 >l,k 3 >l 



N = 4m 2 



K-l 



TZ(K,K,K)+ E a k n KK K(K -k x ,K -k 2 ,K - k 3 



k— (fc^ ,fc2 ) 
k 1 >k 2 >l,k 3 >l 



with 



and 



a k = (_i)fci+fc2+fc3(_ 4m 2)4fc 1+ fc 2+ fc 3 
n r {K — k\,K — k 2 , K — k%) 



NcY 
N F 



ki 



n r (K, K, K) 



(2.20) 

(2.21) 

(2.22) 

(2.23) 
(2.24) 
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for some constants C±, C%, C3. The inequality in (|2.24p follows easily from our counting 
estimates fj2 .111) - fj2. 13f> above. It follows from these expressions that the ratio (12.20P has a 
well-defined limit K — > 00 for small m; indeed, by (|2,18p . one has 



lim \\{K,m) 

K^too 



< (constant) m 2 



(2.25) 



We have considered the maximally tiled graph contributions. It now should not be hard 
to see that counting and adding the contribution of the non-leading non-maximally tiled 
partially reduced graphs in (|2.19|) would only modify (|2.20|) by additional terms carrying 
additional (Nc/Np) suppression factors and higher m powers, and thus not qualitatively 
alter this conclusion. 



It follows from (I2.25P that if we are interested in the small m limit, we need only consider 
the round bracket factor in (|2.19p . i.e., all reduced graphs. In other words, the expression for 
the propagator we need consider is: 



Gk(x, x; m) 



K 

E 



\ 



E J T 

1 a,i c ,i yer r (A,i c ,i) I 



(2.26) 



were now we also include all pure tree graphs, which are reduced graphs by definition. 

For large enough m the sums in (|2.26p are absolutely convergent in the K — > 00 and 
large volume limit. This is the statement of the hopping expansion convergence at large 
m. (The proof, by standard arguments, is actually already implicit in our graph counting 
above.) To continue (|2.26p outside this region requires some resummation procedure. The 
sum (|2.26p can be performed by a self-consistent summation of all reduced graphs (directly in 
the K — > 00 limit) leading to a self-consistent equation for G(x,x) (section [3]) . The solution 
to this equation for large m reproduces, of course, the hopping graph expansion of (|2,26D , 
But the equation obtained by the self-consistent graph summation is meaningful for all m, 
and one may, in particular, examine whether any physical solutions exist as m — > 0. 



3 Self-consistent graph summation and equation for G(x, x) 

All reduced graphs can be generated by iteration of the two basic structures shown in Fig. H] 
(a). Summation over all possible orientations for both the trunks and the attached plaquette 
in Fig. |4] (a) is always implied when attaching these two basic structures at any site. The 
simplest contributions are "one-trunk" graphs obtained by attaching these two basic struc- 
tures at the given site x. Next, there are contributions from n-trunk graphs, any integer n, 
obtained by attaching n of the two basic structures at x. Fig. [5] (b) shows this for n = 2. 
The complete set of reduced graphs can now be generated in a recursive manner by attaching 
these structures again at each available site x + jl, x + ft + i>, x + p, + v + k, x + fi + k on 
top of each trunk of these n-trunk graphs. Indeed, this amounts to starting generating the 
complete propagator G(y, y) at y again, where y stands for any of these sites. Note that the 
subset of pure trees is generated [2] by iteration of the second basic structure of Fig. H] (a) 
above by itself. 
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y + fi + v r, y + ft + O + k 



y + fi + k 



y + jj- 



(a) 



(b) 



Figure 4: (a) The two basic structures attached to a site y whose iteration generates all 
reduced graphs, (b) Attachment of the two basic structures in (a) n = 2 times to the site x 
(two particular directions are shown) to begin generating all "two-trunk" reduced graphs. 



In this recursive rebuilding, however, there can be a slight overcounting of tree segments 
[3], [11] as illustrated in Fig. [5j Though inessential for the conclusion, this overcounting is 
easily corrected by reducing by one (no backtracking) the available lattice directions when 
attaching tree segments at sites beyond the initial point x. This is done by performing 
the recursive rebuilding in two steps [3], |11| : one (i) first recursively builds a complete 
propagator Gi(y,y) at site y incorporating the no-backtracking constraint on tree segments; 
then (ii) obtains the complete propagator G(x, x) by attaching Gj on top of the n-trunk basic 
structures attached at x. 



X + fj, 



■ v 



(a) 



X + V 



X + fl 



X + V 



(b) 



(c) 



y + A 



Figure 5: 2-trunk tree at site x in (b) is contained in the set of 1-trunk trees at x in (a) if all 
2d directions are allowed at x + £i, i.e., backtracking is allowed, (c) Similar ambiguity for the 
attachment of a unit segment at y versus y + (i on top of the first basic structure in Fig. H] 
(a). 



This recursive building of all reduced graphs, subject to the no-backtracking constraint, 
then immediately implies the self-consistency relation for G\ graphically depicted in Fig. [6l 
The n = term stands for the lowest order contribution given by just the bare propagator 
(I2.10p . Note that, with space-time independent source m, Gi(y,y) = Gj, and similarly 
G{y,y) = G, are in fact site-independent by translation invariance. Letting Gj = g(m)l c l F , 
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y 



oo 



n = 




Figure 6: Self-consistent equation for Gi (cf. eq. (I3.27|) ). The hatched blob represent the 
full Gj. 



and using the explicit expressions (12. 6p and (I2.10p the equation in Fig. [6] becomes: 



g[m) 



n.=0 



N F \fg(m)Y (2d-l) f g{m) 



1 

m 



m- 4<f (d- 1) — - 



(3.27) 



N c J \ 2 J 2 V 2 

(|3.27p is to be solved for g(m). The complete propagator G is then given in terms of Gj by: 



G 



m-\Ad\d-l)[jf- 



N F \ fg(m) 



d 



g(m) 



(3.28) 



Had we allowed the slight graph overcounting avoided through (|3.27p - (|3.28[> . we would 
have obtained a self-consistent equation directly for G by now writing G = g(m)l c l F with 
g(m) satisfying: 



g(m) 



m- Ad\d-l) -f- 



N F \ fg(m) 



d 



g(m) 



Clearly, existence of a solution for (|3.27j) implies the same for (|3.28p and vice-versa. 



(3.29) 



4 Absence of chiral symmetry breaking at large Np/Nc 

We consider then solutions to (|3.27p . Letting z = g/2 it assumes the form: 

F(z,m) + 1 = 0, (4.30) 

where 



F(z,m) = 8d 2 (d 



40 



(2d - l)z 2 - 2mz . 



(4.31) 
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For sufficiently large m at fixed (Np/Nc) (|4.30p is immediately seen to have a solution of the 
form 

2z = g(m) = l/m + 0(l/m 3 ) . (4.32) 
This, as noted above, is indeed the expected large m hopping expansion solution. 

We are, however, interested in solutions for small m. Now, F(0, 0) = and F(z, 0) > 
for \z\ > [(2d - l)/8d 2 (d - 1)] 1 / 8 (A^ c /^f) 1/8 ; whereas at its minimum z m = ±[(2d - 
l)N c /md 2 (d - ^Np} 1 / 8 this function has the value 



F(z m ,0) = ~ 
5 



" {2d -I) 5 ' 


" 4 ( 


40d 2 (d-l)_ 





Nf, 



(4.33) 



It follows that F(z, 0) + 1 > for all z, and (|4.30p has no real solutions, provided 



5 



(2d - 1) 
A0d 2 (d- 1) 



1 



For (Nf/Nc) values satisfying (I4.34h then there is no real solution in the absence of sources, 
i.e. m = 0, and no condensate forms. Equality in (|4.34p gives an estimate of the critical 
(N F /N c )c at m = 0. 

Similarly, for m > there is a critical (Np/Nc) c , depending on m, above which there 
is no real solution to (|4.30p . It is immediately seen from (|4.30p - (|4.3ip that this critical 
(Np /Nc) c is an increasing function of m. Thus, for a given (Np/Nc) above the bound (|4.34j) 
there is a mass m for which this (Np/Nc) is critical. A real solution exists for m' > m, which 
approaches (|4.32p for m! large enough, but it turns complex when m' < m. 

In contrast, for (Np/Nc) violating the bound (|4.34|) . (|4.30p has real solutions for all m. 
In particular, in the (Np/Nc) -> limit, (f4"T30]) is solved at m = by g(0) = 2/V2d - 1, 
which, substituted in (|3.28j) . gives 




G =^o) 1 " 1 ' = VdV 1 -^ 101 " (435) 

which in turn, by (|2.3p . yields the well-known result for non- vanishing chiral condensate at 
large (Nc/Np) [2] - [3], [11]. For small m the solution [H] gives a condensate proportional 
to m, whereas at large m it approaches (|4.32|) . 

It is important to note that this m = real solution at small (Np/Nc) does not drop to 
zero as (Np/Nc) is increased to the critical value. Rather, at the critical value, the solution 
turns complex signaling that a condensate ceases to form. This indicates the presence of a 
first order transition, as indeed seen in the Monte-Carlo data [5]. The feature that (Np/Nc) c 
is an increasing function of m is also clearly seen in these data. In this connection also note 
that the method used here, i.e., expansion at large m, resummation and continuation in m, 
can approach this phase transition boundary only from the chirally broken phase. As always 
for two such distinct phases, a different type of expansion would be required to directly 
describe the other, i.e., the chirally symmetric phase. 
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5 Conclusion 



We have considered U(Nq) and SU(Nc) lattice gauge theories at strong coupling ((3 = 0) 
coupled to Np flavors of staggered fermions and investigated the formation of a chiral sym- 
metry breaking condensate within a fermionic hopping expansion. Analysis and classification 
of graphs in their dependence on N@, Np and m led us through a graph resummation proce- 
dure to an equation for the determination of the condensate. At low values of the parameter 
(Np/Nc) this yields the well-known solution for a non- vanishing condensate persisting at 
m = 0. At large values of (Nf/Nc), however, different classes of graphs dominate and the 
situation changes drastically: there is a critical value of (Np/Nc) above which no real solution 
for a condensate exists. The sudden disappearance of a non-vanishing condensate solution 
at a critical (Np/Nc) indicates a first order transition. This result, as well as some other 
qualitative features of behavior as a function of (Np/Nc) an d m are in qualitative agreement 
with the recent simulation results for Nc = 3 in [5]. 

There are some immediate obvious extensions of this work one could consider. Including 
the baryon sector in SU(Nc) in the hopping expansion can be done, but it should not 
materially alter the conclusion. The effect of non- vanishing but small plaquette coupling can 
be incorporated by strong coupling expansion in j3. It is straightforward to incorporate the 
corrections due to this expansion, which is convergent for small enough /?, in our formalism 
above. It is in fact easy to see that doing this has one qualitative effect: it lowers the value of 
the critical (Np/Nc) as a function of increasing (3. It would be interesting to investigate how 
far this can be pushed into the higher (3 region by analogous graph resummations and/or Pade 
approximants, and whether new phase boundaries appear in the (f3, (Np/Nc),m)-space. 

The remarkable feature here is of course the existence of strongly coupled lattice gauge 
systems with colorless ("confined") excitation spectrum and chiral symmetry. This clearly 
opens up several issues, in particular the basic field-theoretic question of the complete phase 
diagram in (f3, (Np/Nc), m ) an d the issue of the existence of continuum limits |15| . There 
are also several areas of potential phenomenological application. One may indeed imagine a 
gauge theory defined at some short scale on a lattice, either physical or serving as a cutoff. In 
strongly correlated electron systems such gauge systems often emerge as effective descriptions 
of underlying electron interactions. Broken chiral symmetry in the gauge system would 
correspond in this case to an anti-ferromagnetic phase. There have been suggestions that 
a chiral symmetry preserving gauge theory, on the other hand, would describe the (high 
frequency regime of) the pseudo-gap phase in the high T c superconductivity phase diagram 
|18j . In particle physics such gauge systems offer a mechanism for obtaining light bound states 
which suggests some obvious potential phenomenological applications to which we hope to 
turn elsewhere. 
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